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We report on experiments exploring the physics of dipolar quantum gases using 
a 52 Cr Bose-Einstein condensate (BEC). By means of a Feshbach resonance, it 
' is possible to reduce the effects of short range interactions and reach a regime 

where the physics is governed by the long-range, anisotropic dipole-dipole in- 
teraction between the large (6 /^b) magnetic moments of Chromium atoms. 
'"^J ■ Several dramatic effects of the dipolar interaction are observed: the usual in- 

' version of ellipticity of the condensate during time-of flight is inhibited, the 

Q ' stability of the dipolar gas depends strongly on the trap geometry, and the ex- 

, plosion following the collapse of an unstable dipolar condensate displays d-wave 

like features. 

Keywords: Bose-Einstein condensation, dipolar quantum gases, Feshbach res- 
^ onances, condensate collapse, vortex rings. 

QO ■ 1- Introduction 



Although quantum gases are very dilute systems, most of their properties 
are governed by atomic interactions. This allows to use them, for example, 



00 

o 

00 ' as quantum simulators to study the many-body physics of systems usually 

encountered in condensed matter physics. 1 However, in all usual quantum 
gases, the interactions can be described extremely well by a short range, 
isotropic contact potential, whose magnitude is proportional to the s-wave 
scattering length a characterizing low energy collisions. 

The dipole-dipole interaction taking place between particles having a 
permanent electric or magnetic dipole moment has radically different prop- 
erties: it is long-range and anisotropic, as one readily sees on the expression 

TT i \ ^ 1 3 cos^ 

^dd(r) = — 73 (1) 

giving the interaction energy £/dd between two polarized dipoles separated 
by r (8 is the angle between r and the direction along which the dipoles are 
pointing). These specific properties have attracted a lot of interest recently, 
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and a large number of theoretical predictions have been made concerning 
dipolar quantum gases (see e.g. Ref. 2 for a review): for instance, the sta- 
bility of a dipolar BEC depends crucially on the trap geometry (see section 
3 below); in a quasi two-dimensional trap, the excitation spectrum can 
display a roton minimum instead of the usual Bogoliubov shape; finally, 
fascinating new quantum phases (including supcrsolids) are predicted to 
occur for dipolar bosons in an optical lattice. 

In practice one always has a competition between contact and dipolar 
interactions; it is therefore useful to define the following (dimensionless) 
ratio of the dipolar and contact coupling constants: 

£dd = 12^V (2) 

The numerical factors are chosen in such a way that a homogeneous dipo- 
lar condensate is unstable against collapse for £dd > 1- For usual atomic 
magnetic moments /i (e.g. for the alkalis), £dd is very small (typically a few 
10~ 3 ) and dipolar effects are extremely small. Here, we report on experi- 
ments with 52 Cr, which has £dd — 0.16 due to its large magnetic moment 
jU = 6^b, and which also allows, via Feshbach tuning of the scattering 
length a, to even enhance £dd- 

This paper is organized as follows. We first describe briefly in section 2 
our experimental setup, with an emphasis on how we use a Feshbach reso- 
nance in order to enhance dipolar effects and create a 'quantum ferrofluid'. 
Section 3 is devoted to the study of the geometry dependence of the stability 
of a dipolar BEC. Finally, we describe in section 4 the dynamics following 
the collapse of an unstable dipolar condensate. 



2. Enhancing dipolar effects using a Feshbach resonance 

A BEC of 52 Cr containing about 50,000 atoms was obtained in 2005 by 
evaporative cooling of optically trapped chromium atoms. 3 Shortly after 
the achievement of condensation, a first effect of the dipole-dipole inter- 
action could be observed in time-of flight experiments: 4 The dipole-dipole 
interaction tends to elongate the BEC along the magnetization direction. 
However, due to the small value of £dd — 0.16, the dipolar interaction was, 
in this experiment, only a small perturbation of the contact interaction, 
which essentially governed the expansion dynamics. 

The existence of several Feshbach resonances 5 in 52 Cr opens the possi- 
bility to tune the scattering length a using an external magnetic field B, 
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Fig. 1. Free expansion of a dipolar condensate for two different orientations of the 
dipoles with respect to the trap axes. The black (resp. gray) squares and lines correspond 
to a situation where the dipoles point along the weak (resp. strong) axis of the trap, 
(a): Perturbative regime e^d = 0.16; (b): £<jd = 0.75. In that case, the dipolc-dipolc 
interaction is strong enough to inhibit the usual inversion of ellipticity in time of flight. 

according to 

Here, ab g — 100a is the s-wave scattering length, B the resonance posi- 
tion, and A the resonance width. The broadest Feshbach resonance in 52 Cr 
is located at Bq = 589 G and has a width A of only 1.5 G. This implies 
that the field control at the level of 3 x 10~ 5 r.m.s. that we implemented 
allows us to tunc a close to with a resolution of about one Bohr radius. 

This 'knob' allowing one to change a allowed us to perform time of flight 
experiments for two different orientations of the dipoles with respect to the 
trap axes, 6 as in Ref. 4, but now for increasing values of £dd- Experimen- 
tal results are shown in Fig. 1. One clearly sees the dramatic effect of an 
increase of £dd on the expansion dynamics. In particular, for £dd — 0.75, 
the inversion of ellipticity of the condensate during time of flight (the usual 
'smoking-gun' evidence for BEC) is inhibited by the strong dipolc-dipolc 
interaction. 

The solid lines in Fig. 1 are theoretical predictions (without any ad- 
justable parameters) based on the Gross-Pitaevskii equation (GPE) gen- 
eralized to take into account the non-local dipole-dipole interaction in the 
description of the macroscopic wavefunction ip(r,t) of the BEC: 

ih^=(-^A + V cxt +g\^\ 2 + jWr\t)\ 2 U dd (r-r')dr^. (3) 
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Here g — 4irh 2 a/m is the contact interaction coupling constant. 

In this set of experiments, the trap geometry was not very far from 
spherical, which limited the study of dipolar condensates to values £dd < 1- 
To go beyond this value and reach the purely dipolar regime e dd 3> 1, 
we shall now see that one needs to tailor the confining potential, so that 
the attractive part of the dipole-dipole interaction does not destabilize the 
condensate. 

3. Geometrical stabilization of a purely dipolar condensate 
3.1. Experimental study 

It is well known that a BEC with attractive contact interactions is unsta- 
ble against long-wavelength fluctuations (this phonon instability leads to a 
collapse of the BEC having a < 0) . As the dipole-dipole interaction has an 
attractive part for dipoles in a 'head-to-tail' configuration — see equation 
(1) for 9 ~ — , it is intuitively clear that in a prolate trap with the dipoles 
pointing along the weak direction of the trap [see Fig. 2(a)], the net effect of 
the dipolar interaction is attractive. Thus, in this configuration, one expects 
that when a is reduced, the condensate becomes unstable, at a critical value 
a cr ;t which should be positive (the small repulsive contact interaction being 
unable, at this point, to counteract the dipolar attraction). Conversely, in 
an oblate trap with the dipoles pointing along the strong confinement direc- 
tion, the critical scattering length should be negative, and a purely dipolar 
quantum gas can be stabilized. 

In Ref. 7, this geometry-dependent stability of a dipolar condensate was 
studied experimentally. A long period (~ 8 um) optical lattice, obtained 
by interfering two laser beams at 1064 nm under a small angle of 8°, was 
superimposed onto the optical dipole trap, allowing us to realize traps with 
cylindrical symmetry around the z-axis (polarization direction) and having 
an aspect ratio A = uj z /uj p that could be varied over two orders of magnitude 
(from ~ 0.1 to ~ 10) while keeping the average trapping frequency CD = 
(uJzUjp) 1 / 3 constant. The experiment consists in creating a BEC in a trap 
with a given aspect ratio A, then ramping a to a final value Of, and finally 
measuring the atom number N in the condensate. One observes that when 
Of is decreased below a critical value a cr it, N suddenly drops to zero. We 
stress that for all the traps we used, the condensate density was roughly 
the same. 

Figure 2(b) shows the measured a cr ;t as a function of A and clearly 
displays the expected behavior: for small A (prolate traps), a cr i t is positive, 
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Fig. 2. (a): Intuitive picture of the geometry-dependent stability of the dipolar Bose 
gas. For a prolate trap (aspect ratio A = lu z /lu p smaller than one) with the dipolcs 
pointing along the weak axis of the trap, the dipole-dipole interaction is essentially 
attractive; such a condensate is thus unstable. For an oblate trap, the dipolc-dipolc 
interaction is essentially repulsive and the BEC is stable, (b): Stability diagram of a 
dipolar BEC in the plane (A, a). The points with error bars are the experimental results 
for the critical scattering length a cr i t below which no condensate is observed; the solid 
line is the stability threshold obtained with a simple gaussian ansatz (see text). 

and starts to decrease when the trap becomes more oblate. For A ~ 10, one 
has a cr ; t ~ 0, meaning that a purely dipolar quantum gas (edd — * oo) can 
be stabilized by an appropriate trap geometry. 



3.2. A simple theoretical model 

A simple way to go beyond the qualitative picture above and obtain an es- 
timate for the instability threshold a crrt (A) is to use a variational method. 
Inserting a Gaussian ansatz (with the axial and radial sizes cr z and a p as 
variational parameters) into the Gross-Pitaevskii energy functional whose 
minimization gives the GPE (3), one obtains the following energy to mini- 
mize: 

NhQ ( 2 1 \ Nhw , , ~ 9N 

e{,„, °.) = — [^ + ^) + ium ( 2 »- + A »•) 

(4) 

where I = yjh/(muj). The first two terms arc the kinetic and potential 
energies, while the third arises from contact and dipolar interactions. The 
function / is monotonically decreasing from 1 to —2 as a result of the 
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anisotropy of the dipolar interaction. For a given A, one can find a (possibly 
local) minimum of E at finite values of (a p ,a z ) if and only if a is larger 
than a critical value: this defines the stability threshold a cr i t (A) within this 
model. 

The solid line in Fig. 2(b) is the result obtained with this simple proce- 
dure, for the experimental parameters lo = 27rx800 Hz and N ~ 2x 10 4 . One 
obtains a relatively good agreement with experimental data. A numerical 
solution of the GPE (3) gives even better agreement with measurements. 8 

Equation (4) allows one to understand easily, in the N — > oo limit, the 
behavior of a cr ;t(A) for A — > and A — > oo. Indeed, in this limit, it is the 
sign of the interaction term which determines the stability; therefore one 
has 

{A — > : BEC unstable if a < add fc .\ 
X — > oo : BEC unstable if a < — 2add, 

where add is the length defined in such a way that £dd = a-dd/a (for 52 Cr, 
one readily calculates, with help of (2), that add — 15a ). It is apparent on 
Fig. 2(b) that for N = 2 x 10 4 , the results are already close to the N — > oo 
limit (5). 



4. <i-wave collapse of a dipolar condensate 

It is natural to ask what happens if one drives the condensate into the 
unstable regime, e.g. by decreasing the scattering length below a cr i t . In 
the case of pure contact interactions, a collapse of the condensate, followed 
by an explosion of a 'remnant' BEC (Bose-Nova), has been observed in 
several systems. 9-12 More recently, the formation of soliton trains has also 
been reported. 13,14 

We have studied the collapse dynamics of a dipolar condensate 15 (in a 
roughly spherical trap) by ramping down rapidly the scattering length to 
a final value of <~ 5a < a cr ;t, then waiting an adjustable holding time, and 
performing a time of flight of 8 ms before imaging the cloud. Figure 3(a) 
presents the evolution of the condensate when the holding time is varied. 
One observes that the cloud, initially elongated along the magnetization 
direction z (horizontal axis on the figure) acquires rapidly a complicated 
structure with a four-fold symmetry, corresponding to a density distribu- 
tion having a torus-like component close to the plane z = and two 'blobs' 
close to the z-axis. Interestingly, this angular symmetry of the cloud is very 
close to the one of a d-wave oc (1 — 3 cos 2 6), i.e. precisely the symmetry of 
the dipole-dipole interaction (1). During the same time period, the atom 
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Fig. 3. (a): Experimental images of the 'exploding' remnant condensate after the col- 
lapse, as a function of the holding time, (b): Result of a numerical simulation of the 
experiment, without any adjustable parameter. The field of view is 130 ^m X 130 fim. 



number in the condensate strongly decreases due to the three-body losses 
occurring because of the high densities transiently reached during the col- 
lapse. 

The group of M. Ueda in Tokyo performed a three dimensional numerical 
simulation of the GPE (3), in which all input parameters were given their 
experimentally measured value. 15 Three-body losses were accounted for by 
adding the imaginary term 



ihL 3 4 



3 body ^ 



M 4 V> 



to Eqn. (3), where L 3 ~ 2 x 1CT 40 m 6 /s is the measured three-body loss 
coefficient. Figure 3(b) represents the results of the simulation. The agree- 
ment is excellent, all the more if one keeps in mind that no adjustable 
parameter is introduced. Let's mention that the inclusion of a small delay 
(also measured independently) in the time variation a(t) of the scattering 
length, due to eddy currents in the vacuum chamber, had to be included 
to achieve a quantitative agreement! The simulation also reproduces quan- 
titatively the time dependence of the condensate atom number. 

A fascinating prediction of the numerical simulation is the spontaneous 
formation, during the collapse, of two quantized vortex rings with opposite 
circulation (and charge ±1), as a result of the strongly anisotropic col- 
lapse: the collapse in the radial directions is fast and quickly followed by an 
outward flow, while axially the flow is still inward, thus giving rise to the 
circulation. Detecting experimentally the presence of vortex rings is very 
challenging, but might be done by using interferometric techniques (e.g. 
matter wave heterodyning) to reveal the winding of the phase of the BEC 
wavefunction around the topological defects. 
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5. Outlook 

The results presented in this paper are the first dramatic manifestations 
of dipolar effects in quantum gases, and pave the way for future studies 
involving even more strongly interacting dipolar systems, especially the 
ones that may be obtained using the permanent electric dipole moments of 
heteronuclear molecules in their ground state. Due to the large value of such 
dipole moments (on the order of one Debye), the long-range character of 
the dipolar interaction could then be used to achieve novel quantum phases 
in optical lattices, 2 as well as to implement promising quantum information 
processing schemes. 16 

However, already in the case of the comparatively weaker magnetic 
dipoles, extremely interesting theoretical proposals deserve experimental 
study; to mention only one example, the generation of two-dimensional 
solitons 17 ' 18 (whose stability arises from the long-range character of the 
dipole-dipole interaction) is a very appealing experiment. 
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